Frontiers in Neuroinformatics

Neurophysiology datasets are becoming increasingly complex, combining behavioral measurements with high-dimensional neuronal activity recordings coming from optical and/or electrophysiological measurements. The Neurodata Without Borders (NWB) standard has emerged in the community as the format of record. While standardized and widely used preprocessing tools generating NWB files have been developed, extensible frameworks for scientific analysis downstream of the NWB ecosystem are still under-represented. We present CICADA, a Python framework dedicated to analysis of neurophysiological data in the standardized NWB format. The toolbox is built as three hierarchically-organized packages: cicada-nwb (NWB access layer), cicada-analysis (plugin-based analysis engine and tool library), and cicada-gui (PyQt5 desktop application at the head of the pipeline). Beyond this architectural separation, CICADA is built around a central design principle: supporting a continuum from turnkey use to full modularity. Researchers can use the complete GUI-driven cicada-gui workflow without writing code, programmatically use existing analysis plugins from cicada-analysis, contribute to new analysis plugins, reuse utilities from cicada-tools, or build entirely custom pipelines on top of the cicada-nwb access layer alone. The same analysis plugin runs identically in interactive GUI and parameter-configured headless modes, enabling reproducible multi-session, multi-animal group analyses. We illustrate the versatility of CICADA with example analyses of behavioral, calcium imaging (two-photon and widefield) and extracellular electrophysiology datasets from rodent laboratories. CICADA is open source, actively maintained, and designed so that any laboratory can contribute at any level of the stack without modifying the core framework.

CT perfusion (CTP) is central to acute-stroke and oncologic imaging, yet quantitative outputs vary substantially across vendor software, undermining reproducibility. We present an open, transparent core (ctp-core) that fits first-pass time-attenuation curves with a gamma-variate model, derives perfusion indices (peak enhancement, time-to-peak, bolus-arrival time, and area under the curve) analytically from the fitted parameters, and renders parametric maps with the ASIST-Japan standardized lookup table (a-LUT) so that visualization is comparable across sites. Every parameter, bound, and processing step is exposed. The method is validated on Monte-Carlo synthetic curves with known ground truth; no confidential or patient data are used. Across signal-to-noise ratio (SNR) levels 5 to 100 (200 independent runs per level) the pipeline recovers peak time to within 0.03-0.52 s and peak amplitude to within 0.4-8.1% (mean absolute error), degrading monotonically with noise; at a representative SNR of 20 it recovers peak time within 0.13 s, peak amplitude within 2.0%, and bolus-arrival time within 0.51 s, with fit quality R-squared = 0.98. The reproducibility demonstration is deterministic (fixed seed) and re-runs to bit-stable metrics. All code, the synthetic-data generator, the standardized-visualization module, evaluation scripts, and a 34-test suite are released openly for independent verification. The contribution is a fully open, parameter-transparent gamma-variate plus standardized-visualization pipeline with reproducible synthetic benchmarks: a reference others can audit, reuse, and build on.

Traumatic brain injury produces widespread axonal damage can be assessed histologically using amyloid precursor protein (APP) immunohistochemistry, which labels injured axonal profiles at cellular resolution [1, 2]. However, quantification of APP pathology remains a major bottleneck: annotation is manual, time-consuming, spatially localized, and variable across raters, limiting scalability and reproducibility. This limitation is particularly important in studies that use histology as a reference for neuroimaging or other tissue-level measurements, where cellular APP pathology must be quantified in a spatial form that can be aligned with imaging abnormalities. Here, we introduce PIGMENT, an annotation-efficient deep-learning framework for automated segmentation and quantification of APP-positive pathology in porcine white matter histology. PIGMENT uses a compact SegFormer-B0 architecture trained on 525 expert-annotated 512 x 512-pixel tiles from four APP-stained sections across three pigs. Because APP-positive profiles are sparse, fragmented, stain-variable, and morphologically diverse, PIGMENT combines limited expert labels with APP-specific augmentation designed to model variation in APP-positive intensity, size, continuity, fragmentation, and local tissue context. We evaluated PIGMENT using an instance-level detection rate that measures whether discrete APP-positive components are localized. Across held-out APP-stained data, PIGMENT achieved a mean instance-level detection rate of 0.86. Across the configurations tested, the highest mean detection rate was achieved by a training set that included sections from different animals, suggesting that annotation diversity may be an important factor under limited-label conditions. By extending limited high-confidence expert annotations into whole-section APP burden maps, PIGMENT provides a scalable framework for characterizing the extent and spatial distribution of traumatic axonal injury. These maps may support future studies that align histological injury burden with imaging-derived measures.

Wave recurrent neural networks (wRNNs) are biologically inspired recurrent architectures that use traveling-wave dynamics to support sequence learning and memory. However, their input-to-hidden pathway remains relatively simple compared with biological neurons, where dendrites perform nonlinear input integration. In this study, we introduce the Dendritic Wave Recurrent Neural Network (DWRNN), which augments the input pathway of the wRNN with nonlinear basal dendritic branches while preserving the original recurrent wave dynamics. We evaluate DW-RNN on a simple copy task, sequential MNIST (sMNIST), permuted sequential MNIST (psMNIST), and noisy sequential CIFAR-10 (nsCIFAR-10). On the copy task, DW-RNN shows learning behavior comparable to the standard wRNN, suggesting that dendritic input integration does not disrupt the recurrent wave-based memory mechanism. On the three sequential image-classification benchmarks, DW-RNN outperforms the standard wRNN, improving accuracy from 97.27 {+/-} 0.15% to 97.82 {+/-} 0.12% on sMNIST, from 96.74 {+/-} 0.17% to 96.92 {+/-} 0.10% on psMNIST, and from 54.30 {+/-} 0.79% to 55.65 {+/-} 0.55% on nsCIFAR-10. In addition to improving mean accuracy, DW-RNN exhibits lower across-seed variability on all three classification benchmarks, suggesting that dendritic input integration may improve the stability of wRNN training. Hidden-activity visualizations further show that DW-RNN preserves the characteristic traveling-wave patterns of the original wRNN. These results suggest that dendritic computation and traveling-wave recurrent dynamics provide complementary mechanisms for biologically inspired sequence learning.

Cortico-thalamo-cortical circuits entail extensive trans-thalamic connectivity between cortical areas, yet their structural organization and function remain poorly understood. Here, the thalamocortical projections of several higher-order thalamic nuclei were characterized by retrograde tracing from two cortical areas, the primary somatosensory (S1) and motor (M1) cortices. Cholera toxin B conjugated with different fluorophores allowed for simultaneous detection of projection neurons targeting S1 and M1. A cell detection pipeline based on neural networks was developed to allow semi-automated analysis of large thalamic imaging volumes to quantitatively infer the spatial distribution of projection neurons in the posterior complex (PO) and the adjacent ethmoid nucleus (Eth), nucleus centrolateralis (CL), nucleus paracentralis (PCN), and the nucleus parafascicularis (PF). The arrangement of neurons projecting to both, primary somatosensory and motor cortices, occurs at different connection strengths and was topographically organized in all nuclei studied. Co-injections into both cortical areas revealed projection neurons with axons branching into both S1 and M1 cortices. Our work introduces a pipeline for semi-automated quantitative analysis of thalamic projection patterns that could be useful for connectivity analyses in general. This approach revealed repetitive anatomical patterns in different thalamic nuclei with regard to projection strength, spatial organization and fraction of projection neurons targeting two cortical areas simultaneously.

The primary motor cortex has a prominent role in transcranial magnetic stimulation (TMS). It is one of few locations that provide directly observable responses, and its physiology serves as model or reference for almost all other TMS targets, e.g., through the motor threshold and spatial targeting relative to its position. It furthermore sets the safety limits for the entire brain. Its easily detectable responses have led to closed-loop methods for a range of aspects, e.g., for automated thresholding, amplitude tracking, and targeting. The high variability of brain stimulation methods would substantially benefit from fast unbiased closed-loop methods. However, the development of more potent methods would early on in the design phase require proper models that allowed tuning and testing with sufficient conditions but without a high number of experiments, which are time-consuming and expensive or even impossible at the needed scale. On the one hand, theoretical researchers without access to experiments miss realistic spatial response models of brain stimulation to develop better methods. On the other hand, subjects should potentially not be exposed to early closed-loop-methods without sufficient prior testing. After all, initial, poorly tuned feed-back as needed for closed-loop operation is known to demonstrate erratic behavior. To bridge this gap, we developed a digital-twin-style population model that generates motor evoked potentials in response to virtual stimuli and includes statistical information on spatial (coil position and orientation) as well as pulse-strength-dependent neural recruitment in the population to represent inter- and intra-individual variability. The population data were extracted from a combination of datasets of healthy and depressed subjects to reach large numbers and wide coverage. The model allows users to simulate different subjects and millions of runs for software-in-the loop testing. The model provides all code open-source to stimulate further development.

Deep brain stimulation has demonstrated its therapeutic potential in modulating pathological oscillations associated with Parkinsons disease and epilepsy. However, its efficacy in treating disrupted theta-gamma phase-amplitude coupling seen in memory-related disorders, such as Alzheimers disease, remains poorly understood. While recent studies have targeted the entorhinal-hippocampal circuit, results remain inconsistent. This discrepancy stems from a lack of mechanistic understanding regarding how stimulation protocols affect this circuit. In this work, we present a reduced multicompartment model of the hippocampal CA1 area that reproduces theta-nested gamma oscillations characteristic of healthy neural activity during memory performance. The model comprises pyramidal, basket and OLM cells with simplified morphologies. We also incorporated CA3-to-CA1 axonal projections, providing a foundational framework for studying how stimulation-induced recruitment of afferent pathways modulates CA1 dynamics. By balancing computational efficiency with anatomical accuracy, our model enables systematic investigation of the effects of electrode placement and orientation, as well as stimulation amplitude and frequency on CA1 neural activity. We demonstrate that the excitatory response in CA1 is primarily driven by the recruitment of Schaffer collateral projections. Overall, this work provides a computationally efficient template for exploring diverse stimulation configurations and could be expanded for developing neuromodulatory strategies to restore physiological network dynamics. Author summaryDeep brain stimulation has shown success in treating Parkinsons disease by suppressing abnormal neural activity responsible for movement disorders. However, when applied to memory-related pathologies, such as Alzheimers disease, the therapeutic outcomes remain unpredictable, ranging from cognitive improvement to impairment. This discrepancy highlights a critical gap in our understanding of how stimulation protocols interact with neural dynamics of the targeted circuits. To address this, we developed a computationally efficient model of the hippocampus, which is involved in memory processes, in order to understand how deep brain stimulation might influence its activity. Our model maintains enough biological accuracy to capture essential memory-related neural activity while remaining lightweight enough for rapid execution and systematic exploration of different protocols. This computational efficiency allowed us to conduct systematic investigations of several stimulation configurations to study their effects on hippocampal dynamics. Overall, this model could provide a useful and computationally cost-efficient tool for exploring the mechanisms of deep brain stimulation and help optimize stimulation protocols aimed at alleviating memory disorders.

We propose a study protocol for routine clinical electroencephalograms (EEGs) from public hospitals, which represents a vast resource for neuroscience research. These non-invasive measures of brain function, paired with rich clinical annotations from large and diverse patient populations, are critical for developing robust artificial intelligence (AI) models and conducting population-level studies. This protocol presents a scalable methodology for curating and harmonizing extensive clinical EEG datasets, encompassing over 40,000 individual studies, to facilitate research applications. Key steps include: (i) integration of raw EEG recordings with corresponding clinical records, including neurological reports, diagnostic codes, and potentially medication data; and (ii) spatial standardization of EEG signals by mapping them to a common brain space defined by functional and anatomical landmarks. The resulting harmonized datasets enable the development of large-scale EEG foundation models, the discovery of novel EEG waveform representations, and the creation of normative "brain charts" for electrophysiological assessment across the lifespan. By enabling standardised, large-scale analyses of real-world clinical EEG data, this protocol supports data-intensive solutions for EEG applications and addresses the challenge of generalising AI models. Our approach promotes the translation of AI tools from research to diverse patient populations, advancing population neuroscience.

Mapping hierarchical brain networks within traditional Euclidean space causes significant structural distortion, undermining neuroimaging diagnostic frameworks. While hyperbolic models like the Poincare ball preserve these nested topologies, they demand heavy computational overhead due to intricate Mobius operations and curved geodesics. This paper introduces a highly efficient non-Euclidean framework for analyzing neurocognitive decline utilizing the Beltrami-Klein ball model. By projecting hyperbolic geodesics as Euclidean straight lines, this approach converts complex distance calculations into simple dot products, radically reducing processing demands. We validated our methodology against state-of-the-art Poincare and Lorentz baselines using datasets for Schizophrenia, Parkinsons Disease, and Alzheimers Disease. The Klein-based framework demonstrates superior performance, delivering both higher diagnostic precision and accelerated processing velocities across all three neurocognitive disorders.

SignificanceDigital Holographic Microscopy (DHM) provides label-free quantitative phase images (QPIs) of living cells and has become a powerful tool for studying cellular morphology and dynamics. While most DHM studies have focused on cell-level analysis, the quantitative characterization of neuronal network organization and maturation from DHM images remains largely unexplored, highlighting the need for dedicated computational approaches. AimWe aimed to develop an automated framework combining deep-learning-based image analysis and graph theory to quantitatively characterize the organization, connectivity, and maturation of neuronal networks in primary rat cortical cultures imaged by DHM. ApproachTwo U-Net convolutional neural networks were trained on manually annotated DHM phase images to segment neuronal cell bodies and neurites. The resulting segmentation maps were used to infer putative morphological connections between neurons and generate graph representations of neuronal networks, referred to as graph fingerprints. A panel of 18 connectomics-inspired graph features was then computed to characterize local and global properties of network organization across four stages of culture maturation. ResultsThe mean area under the receiver operating characteristic curves was 0.98 for cell-body and 0.91 for neurite segmentation, indicating near-perfect identification. Graph-theoretical analysis revealed reproducible topological changes during network maturation in vitro, including increased density, reduced modularity, and progressive network integration. Correlation analysis showed that the 18 graph features grouped into two highly correlated families. A Random Forest classifier identified density and modularity as the most informative descriptors, achieving an accuracy of 87% in classifying maturation stages of neuronal cultures. ConclusionsOur results demonstrate that combining DHM, deep-learning-based segmentation, and graphtheoretical analysis enables quantitative characterization of neuronal network organization and maturation from label-free phase images. This framework provides a foundation for future studies of pharmacological experiments, neuronal network phenotyping, and human induced pluripotent stem cell (hiPSC)-derived neuronal cultures, where quantitative assessment of network organization remains a major challenge.

The corpus callosum is a major white matter bundle responsible for connecting both hemispheres. In mammals, due to a variety of causes, the development of the corpus callosum can be impaired - this brain malformation is known as corpus callosum dysgenesis (CCD). The clinical presentation of CCD varies, with patients exhibiting three morphological phenotypes: agenesis, partial dysgenesis, and hypoplasia. Although the first two presentations are easily detectable on MRI scans, the latter is more challenging, as the structure is fully formed but has a reduced area. In this study, we develop (1) a pipeline to generate synthetic MRI scans with apparent anatomical variation and (2) train a U-Net-based tool to automatically segment the corpus callosum of marmosets in both healthy and disease contexts. Methodologically, a custom script was devised to apply rotation and translation to T1-weighted MRI scans at the volume level. Because the slicing grid remains unchanged, these rigid transformations translate into apparent anatomical variations at the slice level. We compared corpus callosum measurements obtained from automatically segmented masks with those from manually delineated masks. The average Dice score was above 0.90, and the Hausdorff distance was below 0.4 mm. We also stratified our cohort according to phenotype (healthy controls and hypoplastic animals). The magnitude of the effect and the significance level observed between the voxel counts of healthy and hypoplastic animals using manually delineated masks were comparable to those obtained via automatic segmentations. These results show that our pipeline can generate a sufficiently varied training pool to build an accurate U-Net segmentation model with high diagnostic capability.

Compared to traditional gross volumetrics, surface- based models provide greater spatial precision for understanding brain alterations related to developmental, neurological, and psychiatric disorders. Large-scale brain initiatives are combining data from around the world to discover and improve illness- related brain markers. Here, we present a toolkit for 3D brain geometry analysis aimed at addressing key challenges facing large- scale neuroimaging studies. Our framework incorporates scalable methods for multisite data integration, site-specific confound correction, accelerated statistical modeling, interpretable machine learning, and interactive results visualization. The toolkit was tested on data from 21 independently collected study samples participating in the ENIGMA Bipolar Disorder Working Group (N = 3,373). Compared to traditional volume features, we show how subcortical shape measures can be combined across study sites to capture spatially complex differences between diagnostic groups and associations with common treatments. Statistical modeling was accelerated using the Fast and Efficient Mixed- Effects Algorithm (FEMA) and achieved a 16-fold reduction in computation time compared to traditional approaches. Machine learning models showed shape features may provide greater predictive performance over traditional volumes for both diagnostic and treatment prediction tasks, with interpretable weight maps providing insights into the local features driving model performance.

Advances in high-throughput mesoscale microscopy and machine learning-based image analysis pipelines have made unbiased whole-brain imaging widely accessible. However, translating the resulting atlas-mapped datasets into biologically meaningful results remains a substantial barrier owing to their sheer magnitude and complex hierarchical organization. Consequently, reporting structure and analysis methods vary widely across studies, under-mining rigor and reproducibility. To address this, we developed a user-friendly data reduction workflow, HERO (Hierarchy-aware Expression Region Organization), designed to perform hierarchy-aware selection, refinement, ranking, and visualization of whole-brain cell detec-tion datasets. The workflow is customizable to specific needs, requires minimal coding expe-rience, and outputs transparent, curated results. HERO is designed to function as a seamless plug-in within larger-scale whole-brain cell-detection analysis pipelines, providing efficient, unbiased region selection to streamline subsequent statistical analyses and comparative evaluations. Although HERO is developed with mouse cell-detection datasets, it can, in prin-ciple, be applied to any atlas-mapped dataset that contains hierarchical information. In sum, HERO offers a standardized analysis workflow to reduce whole-brain cell-detection datasets, transforming raw regional cell counts into curated results and advancing the effectiveness, interpretability, and accessibility of whole-brain imaging in neuroscience.

Introduction Understanding the mechanisms of human neuromaturation constitutes one of the fundamental questions of neuroscience. While it is well described that large-scale brain maturation is initiated within sensorimotor brain regions and progresses to associative cortex, the underlying developmental neurobiology remains to be fully characterized. Animal models have indicated that cortical inhibitory upregulation might be a driver of neurodevelopment. To investigate the hypothesis that cortical inhibitory upregulation plays a similar role in human neuromaturation, we developed a The Virtual Brain (TVB) based computational model (TVB-Child) to explore potential mechanisms of human neurodevelopment. Material and method We created neurodevelopmental dynamic brain network models capturing neurobiological maturation by using the large-scale brain simulator TVB and fitting brain network models to developmental functional MRI (fMRI) from the Human Connectome Project-Development (HCP-D) data set with 640 subjects with an age range of 6-21 years. Age-dependent trajectories in the fMRI data set were first analyzed by combined group-ICA/Dual Regression extracting subject-specific resting-state networks (RSN). Maturational topographical and topological redistribution of these networks were analyzed by linear and non-linear regression of RSN size and degree and strength centrality. Brain network models were fitted to the fMRI functional connectivity obtained from the HCP-D data set. Hypothesizing that cortical inhibition is a driver of neuromaturation, we analyzed spatiotemporal inhibition parameter gradients in the dynamic brain network model for the hypothesized significant correlations with fMRI RSN maturational trajectories. Results While during development frontoparietal (FP) and default mode network (DMN) grew and exhibited an increase in both degree and strength centrality, becoming dominant network hubs, the attention network underwent network pruning with a decrease in size and node degree. The primary sensory network changed little. For the fitted brain network models, we obtained a high degree of reproduction with correlation coefficients between empirical and simulated functional connectivities ranging between 0.80 and 0.95. Values of the feed forward inhibition model parameter wijFFI representing the strength of regional feedforward inhibitory input exhibited the most significant increase with age within the FP and DMN networks. A less pronounced, but significant, age-dependent increase of the inhibitory parameter values were seen in attention networks and no change within primary sensory networks. Conclusion Our study shows that high order (FP, DMN), attention and primary sensory networks exhibit distinct topographical and topological maturation trajectories. Moreover, brain network modeling revealed RSN-specific age-dependent inhibition trajectories, indicating that the model is able to reproduce and thus support candidate mechanisms of neurodevelopment.

Automated severity rating of rapid eye movement (REM) sleep behavior disorder (RBD) movements would enable multi-night monitoring to detect potentially injurious behaviors and provide objective endpoints for clinical trials. We compared a heuristic classifier using optical flow-derived features against V-JEPA2, a self-supervised video foundation model, for clip-level severity classification (3329 mild versus 284 moderate-to-severe) of in-laboratory video-polysomnography infrared recordings in 86 isolated RBD patients. V-JEPA2 with checkpoint fine-tuning and maximum optical flow-based frame sampling achieved the best performance across both evaluation conditions -- Macro F1 of 0.76 and 93% accuracy in the clip-level split, and 0.68 and 85% in the patient-level split -- outperforming heuristic and domain-specific pretrained models. Clip duration was the dominant heuristic predictor. Whole-night severity scores preserved patient-level ordering despite systematic overestimation, with V-JEPA2 achieving a mean absolute error of 25% versus 52% for the heuristic classifier. These findings establish a foundation for objective, home-deployable monitoring of RBD severity.

Circuit theory has been successfully applied to ecological connectivity modelling, notably via the Circuitscape software, which is typically run locally on a laptop or via a server. For downstream geospatial web applications relying on connectivity analysis, backend infrastructure is required, which can be costly and require advanced data governance. Recent developments in WebAssembly now allow fast C++ or Rust code to be run directly in a sandboxed browser environment for edge computing. We present a WebAssembly/Rust toolset with a geospatial data pipeline and efficient edge-computing implementation of connectivity analysis. This approach may be useful for geospatial modelling software where rasters and memory footprint are small enough for the browser context. Our results show that as expected, Circuitscape solves 1000x1000 raster networks 1-2x faster, but requires further file writes. Accounting for total program runtime, our web implementation can be faster for the given context.

golgi is an open-source platform that takes a peripheral nerve from image to stimulated fiber population through a single graphical interface, with an equivalent scriptable Python API and command-line interface for batch and high-performance use. It integrates promptable image segmentation, automated multi-region tetrahedral meshing, anisotropic finite-element solution of the extracellular field with an explicit perineurium contact impedance, generation of realistic fiber populations and their three-dimensional trajectories, and biophysical activation thresholds through interchangeable backends-- NEURON (via PyFibers) and a GPU-accelerated surrogate (AxonML). Every study exports as an integrity-hashed bundle whose image-to-recruitment provenance is verifiable byte-for-byte. golgi lowers the barrier to in-silico peripheral nerve stimulation modeling for experimentalists and clinicians, using a fully open finite-element stack with no commercial dependencies.

Background: Accurate REM identification is critical for diagnosing REM sleep behaviour disorder (RBD), yet many automated sleep staging systems, especially single-channel EEG models trained on healthy cohorts, do not generalise well to real-life polysomnography (PSG) performed in patients. Objective: To compare a feature-based Random Forest (RF) model tuned for RBD with a state-of-the-art single-EEG deep architecture (AttnSleep), and to assess the impact of cohort adaptation and multimodal inputs (EEG, EOG, EMG, ECG). Methods: Experiments used 89 multi-site in-clinic PSGs (SleepWearables Phase-1) plus 53 MASS healthy controls (mean age 63, std 5 years), with 10-fold cross-validation and out-of-fold evaluation. Model performance was assessed using Cohen's kappa, and attention-based modality analysis was performed to quantify signal contributions. Results: When applied out-of-the-box after training on open-source healthy datasets, both models achieved moderate agreement overall (Cohen's kappa = 0.46), but performance declined in RBD, particularly for REM sleep (AttnSleep Cohen's kappa = 0.19 vs RF Cohen's kappa = 0.44), highlighting limited cross-cohort generalisation. The multimodal model improved overall agreement (Cohen's kappa 0.59 - 0.60) and performance in RBD (Cohen's kappa 0.45 - 0.46), with gains most pronounced in REM (Cohen's kappa 0.45 - 0.49). Attention-based modality analysis identified EEG as the dominant signal, increased EOG contribution during REM, and elevated ECG importance during N3. In RBD subjects, EOG weighting increased relative to non-RBD controls (Delta = +0.081). Guided by these weights, a reduced four-channel EEG model matched full multimodal performance in non-RBD subjects, and adding EOG achieved the best overall configuration (Cohen's kappa = 0.61 overall; Cohen's kappa = 0.48 in RBD) with improved REM classification (53% vs 45% recall). Inclusion of EOG also reduced inter-dataset variability in REM staging. Nonetheless, staging performance in RBD remained lower than in controls, particularly for REM. Conclusions: These results highlight the limited generalisability of minimal-sensor models trained on healthy cohorts, the value of mixed cohort-specific training, and the benefit of multimodal integration and attention-guided channel selection, rather than minimal-sensor approaches alone, for robust clinical sleep staging in pathological populations such as RBD.

Background: Automatic sleep staging algorithms are increasingly applied in clinical and home-based recordings. However, their performance may degrade when transferred to new montages and clinical populations. This is particularly relevant in reduced-channel portable PSG and in disorders such as REM sleep behaviour disorder (RBD), where altered sleep architecture may challenge pretrained models. Objective: To evaluate and compare multiple open-source sleep staging algorithms on a minimal portable PSG setup in controls and patients with and without RBD, and to assess the impact of fine-tuning on clinic-ascertained data. Methods: Six open-source models were applied to 76 subjects recruited from three clinical sleep medicine sites. Performance was assessed using accuracy, F1 scores, and Cohen's kappa, both overall and per sleep stage. Each model was evaluated out-of-the-box and after fine-tuning on clinical data. Results: Out-of-the-box performance varied substantially across models (Cohen's kappa 0.21-0.54). Fine-tuning consistently improved agreement, with the best-performing model (GSSC) reaching Cohen's kappa = 0.58 indicating moderate to good agreement. Performance was highest in controls and lower in patient groups. N3 was the most reliably classified stage across models, whereas N1 remained consistently challenging. REM classification improved after fine-tuning in several architectures but remained model, and subgroup-dependent, particularly in RBD subjects. Conclusion: Fine-tuning substantially mitigates domain shift, updating model parameters to align with new data distributions, when applying automatic sleep staging algorithms to portable clinical recordings. Model architecture influences robustness, with feature-learning approaches demonstrating greater adaptability than fixed-feature models. Despite moderate agreement after adaptation, performance, especially for REM and N1 remains insufficient for fully automated diagnostic use in clinical populations.

Background: Electroencephalographic (EEG) potentials evoked by transcranial magnetic stimulation (TMS) offer a direct window into cortical dynamics. Yet, a systematic exploration of their morphological features, analogous to sensory-evoked potentials, is lacking, especially for stimulation outside the motor cortex. Aim: To obtain region-specific properties of frontal, parietal and occipital networks from the time course of TMS-evoked potentials (TEPs). Materials and Methods: We implemented and applied an automatic procedure to compute peak-to-peak amplitude, peak latency, and inter-peak interval of TEPs recorded from 40 neurotypical subjects stimulated over left occipital (n=25), parietal (n=25), and frontal (n=25) cortices. Results: Occipital TEPs showed the largest peak-to-peak amplitude and longest latency of the first waveform component, independently of stimulation intensity and consistent with the recruitment of a large patch of densely interconnected neurons. Concerning later components, both latency and inter-peak interval systematically decreased along the posterior-to-anterior axis, reflecting progressively faster recurrent dynamics from the alpha-dominated occipital circuitry to the tightly coupled loops between frontal cortex and subcortical structures. Parietal TEPs showed intermediate amplitude and latency measures, consistent with the heterogeneous cytoarchitectonic and connectional organization of the superior parietal cortex. Conclusions: Our findings suggest that TEP morphology is shaped by the distinct properties of the stimulated networks, with early amplitude reflecting the extent of local recruitment and later temporal features tracking the rhythm of recurrent activity. This work offers a mechanistically grounded and practically accessible approach, also released as a Python-based tool, that allows to characterize cortical reactivity across different brain-states and populations.